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tt ■ ABSTRACT 

^; 

I) 

^. 

I 9{, = 6 superconformal Chern-Simons theory was proposed as gauge theory dual to Type IIA 

^ I string theory on AdS4 ® CP^. We study integrability of the theory from conformal dimension 

^ . spectrum of single trace operators at planar limit. At strong 't Hooft coupling, the spectrum is 

O ■ obtained from excitation energy of free superstring on OSp(6|2, 2;R)/SO(3, 1) xSU(3) xU(l) 

• ' supercoset. We recall that the worldsheet theory is integrable classically by utilizing well- 

O ' known results concerning sigma model on symmetric space. With R-symmetry group SU(4), 

00 . 

O ■ we also solve relevant Yang-Baxter equation for a spin chain system associated with the single 

^ ■ trace operators. From the solution, we construct alternating spin chain Hamiltonian involving 

/\ • three-site interactions between 4 and 4. At weak 't Hooft coupling, we study gauge theory 

CZ \ perturbatively, and calculate action of dilatation operator to single trace operators up to two 

loops. To ensure consistency, we computed all relevant Feynman diagrams contributing to the 
dilatation opeator. We find that resulting spin chain Hamiltonian matches with the Hamiltonian 
derived from Yang-Baxter equation. We further study new issues arising from the shortest gauge 
invariant operators TvY^YJ = (15, 1). We observe that 'wrapping interactions' are present, com- 
pute the true spectrum and find that the spectrum agrees with prediction from supersymmetry. 
We also find that scaling dimension computed naively from alternating spin chain Hamiltonian 
coincides with the true spectrum. We solve Bethe ansatz equations for small number of excita- 
tions, and find indications of correlation between excitations of 4's and 4's and of nonexistence 
of mesonic (44) bound-state. 



1 Introduction 

In a recent remarkable development, Aharony, Bergman, Jeffries and Maldacena (ABJM) [[T]| 
made a new addition to the list of microscopic AdS/CFT correspondence [|2j|: three-dimensional 
9l =6 superconformal Chern-Simons theory dual to Type IIA string theory on AdS4 x CP^ [|3J. 
Both sides of the correspondence are characterized by two integer-valued coupling parameters 
N and k. On the superconformal Chern-Simons theory side, they are the rank of product gauge 
group U(A^) X U(A'^) and Chern-Simons levels +k, —k, respectively. On the Type IIA string 
theory side, they are related to spacetime curvature and dilaton gradient or Ramond-Ramond 
flux, all measured in string unit. Much the same way as the counterpart between 9(_ = 4 super 
Yang-Mills theory and Type IIB string theory on AdSs x §^, we can put the new correspondence 
into precision tests in the planar limit: 

N^°°, ^ — > oo with X=— fixed (1.1) 

k 

by interpolating 't Hooft coupling parameter X between superconformal Chern-Simons theory 
regime at X ^ 1 and semiclassical AdS4 x CP"^ string theory regime at X ^ 1. 

In the correspondence between 9l = 4 super Yang-Mills theory and Type IIB string theory 
on AdSs X S^, the integrability structure first discovered by Minahan and Zarembo BH led to 
remarkable progress in diverse fronts of the correspondence [^. It is therefore interesting to 
examine if the new correspondence shows also an integrability structure. The purpose of this 
work is to demonstrate integrability structure inherent to the fA^ = 6 superconformal Chern- 
Simons theory of ABJMB 

AdS/CFT correspondence asserts that gauge invariant, single trace operators in superconfor- 
mal Chern-Simons theory are dual to free string excitation modes in AdS4 x CP'^, valid at weak 
and strong 't Hooft coupling regime, respectively. In particular, conformal dimension of the 
operators should match with excitation energy of the string modes. The 9{^ =6 superconformal 
Chern-Simons theory has SO(6)~SU(4) R-symmetry and contains two sets of bi-fundamental 
scalar fields 7^, 7/ (/ = 1,2,3,4) that transform as 4, 4 under SU(4). Therefore, the single trace 
operators take the form: 

o = Tr(y^iy/^---y^%^)c/j';;.;;^^ 



'Selected but nonexhaustive list of contributions in this subject include ID - lISTI . For a comprehensive mid- 
development review, see lfT3l . 

^Integrability in 5V; < 3 superconformal Chern-Simons theory was investigated by Gaiotto and Yin ll22ll previ- 
ously. Their tentative result indicated otherwise. 



In superconformal Chern-Simons theory, chiral primary operators, corresponding to the choice 
of (|1.2I) with Cj^Zj^ totally symmetric in both sets of indices and traceless, form the light- 
est states. In free string theory on AdS4 x CP^. Kaluza- Klein supergravity modes form the 
lightest states. In this work, we study conformal dimension of single trace operators and iden- 
tify integrability structure organizing the excitation spectrum above the chiral primary or the 
Kaluza- Klein states. 

In section 2, we begin with recapitulating the standard argument for integrability of free 
string on AdS4 x CP^ at ?i ^ oo. Recalling the construction of ll23ll and utilizing the idea 
of [EH, we argue that sigma model on OSp(6|2, 2,M)/[SO(3, 1) xSU(3) xU(l)] supercoset has 
commuting monodromy matrices and infinitely many conserved nonlocal charges. In section 3, 
we begin main part of this work. Guided by earlier development in 5V; = 4 super Yang-Mills 
counterpart, we assume integrability and solve Yang-Baxter equations for R-matrices between 
4 and 4 sites in (11.21) . From corresponding transfer matrices, we then find the Hamiltonian takes 
the form of one-parameter family of 'alternating spin chain', whose variants were studied pre- 
viously in different contexts Il25l - [|29ll . In section 4, we study superconformal Chem-Simons 
theory of ABJM at X —>^ in perturbation theory. Pure Chem-Simons theory is free from any 
ultraviolet divergences since the theory is diffeomorphism invariant and hence topological . 
Once matter is coupled, as in ABJM theory, topological feature is lost and the quantum the- 
ory will receive nontrivial radiative corrections. As such, the single trace operators (11.21) will 
acquire nontrivial anomalous dimensions in general. In three dimensions, logarithmic ultravi- 
olet divergence arises only at even loop orders. Therefore, the first nontrivial correction starts 
at two loops. We compute two-loop operator mixing and anomalous dimension matrix of the 
single trace operators (|1.2I) . In dimensional reduction method, we compute the complete set of 
relevant Feynman diagrams and find that the two-loop anomalous dimension matrix matches 
with the integrable 'alternating spin chain' Hamiltonian derived in section 3. In section 5, we 
study a new important feature of the superconformal Chem-Simons theory compared to fA^ = 4 
super Yang-Mills theory. Since the anomalous dimensions begin to arise from two loops and 
next-to-nearest sites, the shortest single trace operators of L = 1 will be subject to 'wrapping 
interactions'. The 'altemating spin chain' Hamiltonian does not describe spectmm of L = 1 
operators, so we compute all relevant 'wrapping interaction' diagrams and construct the cor- 
rect Hamiltonian for L = 1 . Curiously, we find that the correct spectrum coincides with the 
naive spectrum computed from the 'altemating spin chain' Hamiltonian at L = 1. In section 6, 
utilizing results previously obtained for general A„_i Lie algebras [[30ll26ll29ll . we explicitly 
write down eigenvalues of the transfer matrices and Bethe ansatz equations of the 'alternating 
spin chain' we derived in section 3. To gain understanding how the 'altemating spins' behave, 
we solve the equations for a few simple situations. We find an indication for real- space cor- 



relations between excitations on 4 spin sites and those on 4 spin sites, and for non-existence 
of meson-like (44) bound-states. We discuss various implications of these findings for general 
excitations. In particular, we argue that general excitations are more complex than the pattern 
emerging from closed SU(2) sub-sectors discussed recently lf3Tll32ll . 

• Note added : While bulk of this work was completed, we received the preprint by Minahan 
and Zarembo |[33l . which deals with issues overlapping with ours. Version 1 of their preprint 
was based on ungrounded Feynman rules and incomplete set of contributing Feynman diagrams. 
In light of utmost importance of precise evaluation, we decided to carefully carry out various 
internal consistency checks before releasing our results. Meanwhile, we also received preprints 
by Arutyunov and Frolov [|34ll . which overlaps with section 2 and substantiate several pertinent 
issues including kappa symmetry. We also note the preprints by Stefanski ll35ll and by Gromov 
and Vieira [[36l on the same issue. 

2 Integrable String from Worldsheet Sigma Model 

In this section, we set out a motivation for searching for integrability m9\C —6 superconformal 
Chern-Simons theory. The ?i — > oo dual of this theory is Type IIA string theory on AdS4 (8) CP^. 
The background is a direct product of symmetric spaces, AdS4 and CP"^. It is well known that 
the (1 + 1) -dimensional sigma model on symmetric space is classically integrable. So, at least 
for bosonic modes, we expect worldsheet dynamics of a free string on AdS4 (8) CP^ is integrable 
at the classical level, ?i — > oo. In this section, we recapitulate this argument for the bosonic part 
and discuss how the construction to full superstring can be made. 

Bosonic part of string worldsheet Lagrangian on AdS4 ® C¥^ is given by 

(2.1) 



4 = ?z ^ V^h'^^ [{DaXnHDpX"') + (DaZ«)t(DpZ«) 



4%. 



Here, we use embedding coordinates in M?'^ and C^ and describe AdS4 =SO(3,2)/SO(3,l) and 
CP3 =SU(4)/SU(3)xU(l) as G/H coset hypersurfaces: 

AdS4 : (X'") = {X-\X\X^,X^,X^) with X^ = 1 , 

CP^: (Z'') = (Z1,Z2,Z3,Z4)/{~,C} with |Z|2 = 1, (2.2) 

respectively. The hypersurface conditions are imposed by introducing auxiliary connection 
KaAa and by defining covariant derivatives B DaX"" = daX"" + iKaX"" and DaZ" = daZ" + 
iAaZ". These conditions imply that X'"daX„, = and (DaZ")'''Z" = Z"^(DaZ") = 0. Fol- 
lowing [|23ll . we first recapitulate basic aspects for classical integrability of sigma model on 



We introduced auxiliary connection Ka to treat AdS4 in complete parallel to i 



AdS4 X CP^. Construction of the coset sigma model is facilitated by the coset elements: 

G{a) = g{a) ®g{o) = e''^^^^) © e''^^^'') (2.3) 

where P{<5),P{o) are projection matrices onto respective one-dimensional subspaces. They are 

P'""{a) = X"\o)X"{o) with 8,nnP'""{o) = l, 

p^\o) = Z"'^{(5)z''{o) with dabP''\o) = l, (2.4) 

respectively. By elementary algebra, we verify that 

G{g) = G-\(5) = (l5-2P(o))©(l4-2P(a)). (2.5) 

Then, because -S\DaX"'\^ = Tr{dag ■ dag^^) for AdS4 and +4\DaZ"\'^ = M^ag- ^^g^^) for 
CP^, the worldsheet action (12.11) is expressible as 



/bosonic = Y y v/^/?"P [ - -Triaip + 2Tr7a7p 



(2.6) 



where J = g ^dg and J = g Mf, respectively. We shall choose the conformal gauge \J —hh^^ = 
5"P on the worldsheet. This leads to Virasoro gauge condition 

r± = -^{Jo±Ji)^ + iJo±Jif = . (2.7) 

The currents 7, / are conserved by equations of motion, and define tangent flows on the G/H 
coset space. 

We now take group conjugation and transform the left-invariant currents 7, 7 to the right- 
invariant currents: {j,j) = {g-J-g^\ g-J - g^)- The equations of motion in conformal gauge 
are 

d*;(o)=0 and d* J(o) = 0. (2.8) 

From the Bianchi identities, we also have 

dj + jAj = and jJ+JaJ^O. (2.9) 

Finally, Virasoro constraints are 

-^(jo±Ji)' + (/o±Ji)' = 0. (2.10) 

We can solve these equations using the Lax representation. Consider the Lax derivative with 
flat connection a{x) depending on a spectral parameters: 

D{x)=d + a{x) with da + aAa = 0. (2.11) 



Using (|2.8[|2.9I ). we find that the most general form of the Lax connection is given by 



2 2x 



flW = ^-rX<') + ^-r*X<') xeC + W/{±i}. (2.12) 



and similarly construct a{x) from i{o). With the flat connection A (x) = {a{x),a{x)), consider 
the Wilson line 

W[r,x] = 2'exp r Ja{x)\ (2.13) 

As the connection A{x) is flat, the eigenvalues of the Wilson line are independent of the choice 
of the contour y. Thus, all Wilson lines commute each another and provides classical R-matrices 
obeying Yang-Baxter equations. Expanding in spectral parameter x, we then obtain infinitely 
many conserved nonlocal charges as moment of the power series: 



n\ Jy 



Y 

3 



(2.14) 

x=0 



This establishes that the sigma model on AdS4 x CP is classically integrable. 

We now discuss how the above consideration may be extended to Type IIA string on the 
supercoset: 

G ^ OSp(6|2,2,R) 

H SO(3,l)xSU(3)xU(l)' ^ ■ ^ 

With the coefficients of current bilinears determined as in (|2.6I) . we see immediately that the 
worldsheet Lagrangian is expressible as supertrace over the supergroup OSp(6|2,2,M): 

i?2 



/supercoset = y y Str (7 A V) . (2.16) 

Here, 7(a) = G^^(a)dG(a) and G{(5) = exp(/7rP(o)) is the supercoset element. This indicates 
that the bosonic action (|2.6I) is extendible straightforwardly to a supercoset action by adding 
24 fermionic off-diagonal components to (I2.3H2.5I) and define super-projection matrix P and 
supercoset element G analogously. 

Construction of infinitely many nonlocal currents requires a new condition to the supercoset. 
If the supergroup G permits Z4 grading under which the subgroup 7/ is a fixed point set, the 
construction of [i24il implies that a flat connection exists from which nonlocal currents can be 
constructed through the Lax formulation. From the embedding we constructed of, we have 7 = 
J + Q, where Q denotes fermionic current. For the supergroup we deal with, G = OSp(6|2, 2), 
it is well known that G admits no outer automorphism of order four ll37l . However, one can 
easily construct a suitable Z4 inner automorphism. Since we need the subgroup 7/ is a fixed 
point set, the automorphism can be defined as a product of two Z2 involutions on the defining 



representations of SU(4)~SO(6) and Sp(4) modulo overall reflection. This then ensures that 
the G/H current J = Qi ® J ® Q^, is Z4 graded as [1,2,3] and that infinitely many conserved 
nonlocal currents can be constructed accordingly. 

At quantum level, the supergroup G= OSp(6|2,2) has another nice feature that its Killing 
form vanishes identically. This means that sigma model on G would be conformally invariant, 
at least, at one loop. We actually need to quotient G by bosonic subgroup H and consider string 
worldsheet action on the supercoset G/H. This action in general breaks the conformal invari- 
ance. To restore the conformal invariance, a suitable Wess-Zumino term needs to be added. It 
was observed [[38l that the requisite Wess-Zumino term can be constructed provided the bosonic 
subgroup // is a fixed point set of the Z4 grading of G. This is precisely the same condition that 
ensures the existence of a flat connection and infinitely many conserved charges thereof. There- 
fore, the supercoset sigma model is conformally invariant and can describe consistent string 
worldsheet dynamics, at least at one loop order in worldsheet perturbation theory. 

Given such mounting evidences, it is highly likely that Type IIA string on AdS4 x CP is 
integrable at ?i ^ 00 and further extends to X finite and even to weak coupling regime Q With 
such motivation, we now turn to the main part of this work and investigate integrability at the 
weak coupling regime, ?i ^ 0. 

3 Integrable Spin Chain from Yang-Baxter 



The U(N)xU(N) invariant, single-trace operators under consideration 

0«(,)^Tr(y^iyt...y4yt) 

- 0(y)(') = Tr(y;/^i---y,V^) (3.1) 

are organized according to SUie(4) irreducible representations. Operator mixing under renor- 
malization and their evolution in perturbation theory can be described by a spin chain of total 
length 2L. What kind of spin chain system do we expect? In this section, viewing the operators 
(13.11) as a spin chain system and utilizing quantum inverse scattering method, we shall derive 
spin chain Hamiltonian. 

As is evident from the structure of operators (|3.1I) . the prospective spin chain involves two 
types of SU^(4) spins: 4 at odd lattice sites and 4 at even lattice sites. It is thus natural to expect 
that the prospective spin chain is an 'alternating SU(4) spin chain' consisting of interlaced 4 
and 4. To identify the spin system and extract its Hamiltonian, it is imperative to solve inho- 
mogeneous Yang-Baxter equations of SU/;(4) 91-matrices with varying representations on each 



■^We note that the no-go theorem of Goldschmidt and Witten 11391 for quantum conservation laws is evaded in 
the present case since the isotropy subgroup is not simple, and may lead to quantum anomalies ll40l . 



site. In fact, a general procedure for solving Yang-Baxter equations in this sort of situations is 
already set out in ll25l . By construction, resulting spin chain system will be integrable. In this 
section, we shall follow this procedure and find that the putative SU(4) spin chain is an 'alter- 
nating spin chain' involving next-to-nearest neighbor interactions nested with nearest neighbor 
interactions q 

We first introduce y\^{u) and y{^{u), where the upper indices denote SU(4) representations 
of two spins involved in 'scattering process' and m,v denote spectral parameters. We demand 
these R-matrices to satisfy two sets of Yang-Baxter equations: 



^1t{u-v)m^,{u)m^,{v)=m^,{v)Rli{u)Rtt{u-v) (3.2) 

^tiiu-v)D\Uiu)D\f,{v)=D\f,{v)DKtiiu)9Ktiiu-v) (3.3) 






Here, the lower indices z, j denote that the 9^ matrix is acting on i-th and 7-th site V, Vj of the 
full tensor product Hilbert space Vi ® V2 ® ■ • • ® V2L- We easily find that the R-matrices solving 
(13.21 13.31) are given by 



}-44/,,\ ,,TT I 11 „„^ co-44 



Dr*{u) = ul + F and DT'^iu) ^ -{u + 2 + a)I + K (3.4) 

for an arbitrary constant a. Here, we have introduced identity operator I, trace operator K, and 
permutation operator P: 

ihiYiX = s^sj {^mYiX = s^^s^y. imXi = ^'lA , (3.5) 

acting as braiding operations mapping tensor product vector space Vk Vf to itself. 

We also need to construct another set of R-matrices ^{'^^{u) and D\'*'*{u) generating another 
alternative spin chain system. We again require them to fulfill the respective Yang-Baxter equa- 
tions: 

^i{u-v)9lU{u)^§{v)=mf,ivMi{uMi{u~v) (3.6) 

^t^(i.-v)lHf^(z.)<(v) =inl^(v)lHt|(z.)9^t^(..-v) (3.7) 

Again, we find that the R-matrices that solve (13. 6[ 13.71) are given by 



m^{u) = uI + ¥ and 9\^{u) = -{u + 2 + a)I + K, 



(3.8) 



where a is an arbitrary constant. 

^For a construction in SU(3), see ll28l . Generalizations to arbitrary Lie (super)algebras and quantum deforma- 
tions thereof were studied in 



In the two sets of Yang-Baxter equations, the constants a, a are undetermined. We shall 
now restrict them by requiring unitarity. The unitarity of the combined spin chain system sets 
the following conditions: 

m'^{u)m'^{-u) =p{u)i 

m^'\u)m^{-u) =p{u)i 

m'^{u)m^{-u) =o{u)i (3.9) 

where p{u) = p(— m),p(m) = p(— m),o(m) are c-number functions. It is simple to show that the 
first two unitarity conditions are indeed satisfied for any a, a. It is equally simple to show that 
the last unitarity condition is is satisfied only if a = —a. Without loss of generality, in what 
follows, we shall set a = —a = 0. 

Viewing (13.11) as 2L sites in a row, we introduce one transfer T-matrix 

ro(M,a)=in^t(M)9^^^(M + a)9^^^(M)9lS(" + «)---9^02L-i(")^02L(" + «), (3-10) 

for one alternate chain and the other T-matrix 

ro(M,a)=9^^t(M + fl)9ig(M)9^J^(M + a)9^g(M)---lK^^^_i(M + a)9^g^(M), (3.11) 

for the other alternate chain, where we introduce an auxiliary zeroth space. By the standard 
'train' argument, one can show that the transfer matrices fulfill the Yang-Baxter equations, 

^QQ,{u-v)To{u,a)T(y{v,a) = Tff{v,a)To{u,a)^QQ,{u-v) , (3.12) 

and 

D\QQi{u — v)To{u,d)TQf{v,d) = TQi{v,d)To{u,d)9\QQi{u — v). (3.13) 

In addition, by a similar argument, one may verify that 

9^qq,(m — v + a)ro(M,a)ro'(v, — a) = T(y{v,—a)To{u,a)y{QQ,{u — v + a) . (3.14) 

We also define the trace of the T matrix by 

x''^\u,a)=TrTo{u,a). (3.15) 

and 

x^^'(M,a) =TrTo(M,a) (3.16) 

where the trace is taken over an auxiliary zeroth space. 



It then follows from the Yang-Baxter equations that 

[x"'^(a,a),x^'^(v,a)]=0 

[r\u,a),r\v,a)]=0, (3.17) 

and 

[x"i'(M,a),x"i'(v,-a)]=0. (3.18) 

Here, in the first two equations, a, a are arbitrary and denote two undetermined spectral param- 
eters. These parameters are restricted further if we demand the last equation to hold. Indeed, 
the two alternating transfer matrices commute each other if and only if a = —a. 

As for all other conserved charges, the Hamiltonian is obtained [jby evolving the transfer 
T-matrix infinitesimally in spectral parameter u: H = dlogx(M,a)|;,=o where d = d/du. By a 
straightforward computation, we obtain the 44 spin chain Hamiltonian as 

Hie-x = -{2-a)l-{A-a^)¥2i-i,2t+\ 

- (a-2)P2^-l,2^+llK2f-l,2f+(« + 2)P2£-l,2£+lIK212£+l, (3.19) 

where we scaled the Hamiltonian by multiplying (a^ — 4). 

By the same procedure, we also find that the Hamiltonian for for the 44 spin chain is given 
by 

H2i = -(2 + a)I-(4-a2)P2£,2m 

+ (a + 2)P2i',2f+2lK212£+l -{a- 2)P2^,2£+2lK2f+l,2^+2 , (3.20) 

where we have replaced a by a using the relation a = —a. 

At this stage, any choice of the parameter a is possible in so far as hermiticity of the Hamil- 
tonian is satisfied. The latter condition requires that a is a pure imaginary number. Physically, 
we are interested in the situation where 4 ^^ 4 is a symmetry. This is nothing but requiring 
charge conjugation symmetry, equivalently, reflection symmetry in dual lattice. We thus put 



''The following derivation of Hamiltonian is valid only for L>2. This means that the energy eigenvalues of 
the following Hamiltonians for the case L = 1 do not agree with true energy eigenvalues. 



a = zOlI Adding the two alternate Hamiltonians, we get total Hamiltonian|3 : 

2L 

^total = 22^^/+ 1,^+2 (3.22) 

£=1 



with 



H> 



£,£+1^+2 



41 - A¥u.+2 + 2P^,^+2K£ .+1 + 2¥ij+2^i+u 



(3.23) 



In this derivation, there is always a freedom of shifting ground state energy by an arbitrary 
constant. From the outset, we assumed integrability but, except that the symmetry algebra is 
SU/e(4) and that spins are 4,4 at alternating lattice sites, we did not utilize any inputs from 
underlying supersymmetry. With extra input that that supersymmetric ground-state has zero 
energy, one can always fix the freedom. The (|3.23l) is the Hamiltonian after being shifted by +6 
per site accordingly. 

4 Integrable Spin Chain from Chern-Simons 

In this section, we approach integrability from weak 't Hooft coupling regime of the supercon- 
formal Chem-Simons theory. We use perturbation theory and look for a spin chain Hamiltonian 
as a quantum part of the dilatation operator acting on the single trace operators. As mentioned 
above, in three-dimensional spacetime, general power-counting indicates that logarithmic di- 
vergence arises only at even loop orders. Therefore, leading-order contribution to anomalous 
dimension starts at two loops. In general, as well understood from general considerations of the 
renormalization theory, the divergence in one-particle irreducible diagrams with one insertion 
of a composite operator contain divergences that are proportional to other composite operators. 
Therefore, at each order in perturbation theory, all composite operators whose divergences are 
intertwined must be renormalized simultaneously. In addition, renormalization of elementary 
fields needs to be taken into account. This leads to the general structure: 

0t,{Y^,,,jlj = Y^Z''j,0^,,{ZY,,,,ZY,l) (4.1) 

A' 



^Alternatively, one may relax hermiticity of the Hamiltonian and only demand symmetry under parity and time- 
reversal, leading to so-called PT-symmetric system pT]. This again sets a to zero. Strictly speaking, however, this 
latter condition is weaker than the hermiticity requirement. 

^We remark the following useful identities 

Vi,e+2^i.e+i ^^M,e+2Vi,e+2, Ve,e+2Ke+i/+2 = ^e,i+iVej+2 ■ (3.21) 

We shall find them useful later when investigating issues concerning wrapping interactions. 
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For the operators we are interested in, this takes the form of 

<e-IZ^^(A)0l (4.2) 

N 

with the UV cut-off scale A. Therefore, the anomalous dimension matrix A is given by 

dlogA 

hi the rest of this section, we compute anomalous dimension matrix for the single trace 
operators that were associated with the 'alternating spin chain' in the last section: 



O 



v) T'^ / v^i vt vhv't v^Lvt 



(i) - ^M ^ Vi^ ^2 ^ Vz, 



Tr ( y^i YJ Y'^ yJ ■■■ Y^^YJ ) . (4.4) 



In iV; = 6 superconformal Chern-Simons theory, the scalar fields 7^, Yj are bifundamental fields 
of U(N)xU(N) gauge group, and transform as 4 and 4 of SU/?(4) R-symmetry group. In Ap- 
pendix A, we explain field contents and action of the theory in detail Q Schematically, the 
action of the ABJM theory takes the form 

/ = / A (CS{A) - CS(A)) - Tr{DY)]DY^ + Tr ^^^/^^y -V^-Vb. (4.5) 

Here, the Chem-Simons density is given by 

CS(A)=£'""PTr 



2/ 



(4.6) 



Covariant derivatives are denoted as D^, while self-interactions involving bosons and fermion 
pairs are denoted by Vb, V'f, respectively. See Appendix A for their explicit form. We will recall 
them at relevant points in foregoing discussions. 

To extract the dilatation operator, we compute the correlation functions 

/o{jJTr(y^'[y^i---y/^y^^)\ for l^oo (4.7) 

by summing over all planar diagrams in perturbation theory in 't Hooft coupling X. 

In evaluating so, there arises an important issue regarding consistency of regularization with 
gauge invariance and 9\i = 6 supersymmetry. We shall adopt dimensional reduction method 
(See, for example, discussions in [|43l ). This method retains £"'"^ and Dirac matrices always 
three-dimensional. In each Feynman integral, we then manipulate the integrand until all e"'"^ 
and Dirac matrices are eliminated and the integral is reduced to a Lorentz scalar expression. We 
then employ dimensional regularization and evaluate the integral. Still, this leaves out infrared 
divergences that would have been absent were if the theory four-dimensional. As we will be 



^We closely follow notation and convention of f^2\ . 
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only concerned with logarithmic ultraviolet divergences, we will take a practical approach that 
we regularize infrared divergences by introducing mass terms in evaluating Feynman integrals 
in dimensional regularization. We then remove the regulator mass first and then take the space- 
time dimension to three. Previously, it was checked that the dimensional reduction method is 
consistent with Slavnov-Taylor-Ward identities. Yet, to date, it is not known if the method is 
compatible with fA^ = 6 supersymmetry. Thus, in our computations, we shall not assume a 
priori any input related to supersymmetry. Rather, we will put our result to a test against vari- 
ous consequences of supersymmetry — for instance, vanishing anomalous dimensions of chiral 
primary operators and superconformal nonrenormalization theorems. 

Using the convention and Feynman rules explained in appendix, we computed all two-loop 
diagrams that contribute to anomalous dimensions of elementary fields Y',Yj and composite 
operators 0^,1. Acting on the space of the operators, each Feynman diagram can be attributed 
to the braiding operations I, K, P introduced in (13.51) and their combinations. At two loops, we 
computed the complete set of Feynman diagrams that contribute to each of these operators. The 
result turned out 



2L I I 



=1 



2 '''^" ''''^' 2 



(4.8) 



and this is precisely ^ times the alternating spin chain Hamiltonian (13.231) we derived from 
SU(4) Yang-Baxter equations in the last section. In the rest of this section, we explain essen- 
tial steps for deriving the Hamiltonian and relegate technical details of evaluating Feynman 
diagrams in the Appendix. We find it convenient to organize contributing Feynman diagrams 
according to the number of sites that participate in the Hamiltonian. 

• Three-site scalar interactions: 
A salient feature of the alternating spin chain Hamiltonian we extracted in section 3 from cou- 
pled Yang-Baxter equations is that it contains interactions up to next-nearest-neighbor sites. We 
thus need to see if such interaction arises from superconformal Chern-Simons planar diagrams 
and, if so, if the interactions are of the same type. From the Feynman rules (see Appendix A), 
it is evident that scalar interaction — Vb in (14.51) is the source of three-site interactions, whose 
explicit form is given by 

2 



1 /27r , 
Vb = -3 ylTr 



y^y-'yJy^y^y' + Yj'y'yJy-'y^y^ 



+AY/Y-^YIY^YJY^ - 6Y/ Y'yJ Y^^Y^Y 



(4.9) 



The two-loop Feynman diagram is depicted in Fig[T] From planar diagram combinatorics of 
gauge invariant operators at infinite length 2L -^ oo, we find the following contributions arising: 

12 




Figure 1 : Two loop contribution of scalar sextet interaction to anomalous dimension of . 




k + l O 




I . 

li. y k + l O' 




Figure 2: Two loop contribution of gauge and fermion exchange interaction to anomalous di- 
mension of . 

^ll+i +'^£+11+2 from the first two terms, Fg^£j^^2 from the third term, and I-\-F£^gj^^2^e,£+i + 
^£,£+2^1+ 1, £+2 from the last term. Taking account of combinatorial multiplicities, we find that 
the scalar sextet potential contributes to the dilatation Hamiltonian as 



2L 



^B = ^' E ^I - ^£,£+2 + ^^£,£+2^£,£+l + \^£ 



i^e- 



'+1, 



2^^/+i 



(4.10) 



(see Appendix B2). Evidently, compared to the anticipated alternating spin chain Hamiltonian, 
we have discrepancy in on-site (proportional to I) and nearest neighbor (proportional to K) 
terms. These are interactions that would arise from gauge or fermion-pair exchange interactions 
and from wave function renormalization of elementary fields Y,YK 

• Two-site gauge and fermion interactions: 
The scalar fields Y^,Yj are bifundamentals of U(N)xU(N). Their gauge interactions can be 
read off from covariant derivatives: 



L^mJ^ — 0,7j/ + lAfjiY li Af, 



and 



't , .-T 



D„,Y] = d,nY] + iA,nY]-iY]A„, . (4.11) 



As usual, there are paramagnetic interactions (minimal coupling) and diamagnetic interactions 
(seagull coupling). We see that gauge interactions contribute to two-site terms for both I and K. 
Two relevant Feynman diagrams are (a) and (c) in Fig.[2l 

The Feynman diagram contributing to I operator arises from square of diamagnetic interac- 
tions in ^-channel. See Fig. [2la). This diagram is infrared divergent for each subgraphs. We 
regulate them by giving a mass to internal propagators. Upon removing the regulator mass to 
zero, we find a finite part. However, this part turned out ultraviolet convergent and hence does 
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not contribute to anomalous dimension. The Feynman diagram contributing to K operator arises 
from product of diamagnetic interaction and two paramagnetic interactions. See Fig. [2c). Tak- 
ing the net momentum of to zero, which is sufficient for extracting anomalous dimension, we 
find that only one orientation of diamagnetic interaction vertex yields nonvanishing result. For 
details of Feynman rules of gauge interactions and Feynman diagram evaluation, see Appendix 
B3. We found that gauge interactions contribute to the dilatation operator by 

2L 



"gauge — ^ 7 . 



il-^Ko+i 



(4.12) 



Consider next two-site terms induced by fermion-pair exchange diagrams. The relevant part 



of the Lagrangian in (14.51) is the fermion-pair potential: 



271/— 

Vf = — Tr 
k 

Ini 



y^ty/vj/Vvj/^ - lY^Y^m^'^'j + z"^^Y]mjYl^'L 



y/y^txj,^Xj,V - IY^Y]^!^^^ + ZuklY^^'^^Y'^^'^^ 



(4.13) 



From Feynman rules, we see that planar diagrams formed by square of the second terms in both 
lines in (14.131) give rise to K interactions to the two-loop dilatation operator. See Fig. |2b) for 
the relevant Feynman diagram and Appendix B3 for the details of computation. 

In fact, at planar approximation, there is no other Feynman diagrams that contribute to two- 
site interactions [^j. Taking account of numerical weights in (14.131) . we find that the fermion 
potential contributes to the dilatation Hamiltonian as 

2L 

//F = ?t'£K«+i. (4.14) 

1=1 

• One-site interactions: wave function renormalization 

Adding up all the two- site interactions to the three- site interaction, we see that terms involving 
K operator cancel out one another. On the other hand, terms involving I operator add up to 
(1/4)X^. So, up to overall (volume-dependent) shift of the ground state energy, the dilatation 
operator agrees with the alternating spin chain Hamiltonian we derived in the previous section. 
As we are dealing with superconformal field theory, spectrum of dilatation generator bears an 
absolute meaning. Moreover, there could be potential clash between dimensional reduction 
we used and superconformal invariance. Therefore, to ensure internal consistency of quantum 
theory, we shall now compute terms arising from wave function renormalization of 7, 7^. These 
are all the remaining contributions to anomalous dimension of composite operator . 

Wave function renormalization to 7,7^ arises from all three types of interactions. Even 
though there are huge numbers of planar Feynman diagrams that could potentially contribute 



'"For L = 1, however, there will be wrapping interactions. We will discuss them in detail in the next section. 
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(b) 




Figure 3: Two loop contribution of diamagnetic gauge interactions to wave function renormal- 
ization ofY.Y^ . They contribute to I operator in the dilatation operator 





Figure 4: Two loop contribution of paramagnetic gauge interactions to wave function renormal- 
ization ofY.Y^ They contribute to I operator in the dilatation operator 



to wave function renormalization, a vast number of them vanishes identically or cancel one an- 
other. First, diagrams involving gauge boson loops either vanish because of parity-odd nature of 
the gauge boson propagators or cancel among U(N) and U(N) diagrams [Hj. Nonzero contribu- 
tion arise only from diamagnetic interactions shown in Fig. [3l from paramagnetic interactions 
shown in Fig. HI and from Chem-Simons cubic interactions shown in Fig. [51 

Second, diagrams involving vertices in the first and the second lines in Vf (|4.13l) cancel 
by combinatorics and relative coefficients. Hence, the cancellation is attributable to fA^^ = 6 
supersymmetry. The only surviving diagram arise from cross term of vertices in the last line in 
(14.131) . The Feynman diagram is shown in Fig. [6l 

Third, there are also contributions coming from gauge-matter interactions. Again, almost all 
diagrams vanish because of parity-odd nature of gauge boson propagator. The only surviving 



'^ Notice that gauge boson propagator for U(N) and U(N) gauge groups have weight +k and —k, respectively. 




Figure 5: Two loop contribution of Chem-Simons interaction to wave function renormalization 
ofY, Y\ They contribute to I operators in the dilatation operator 
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k + l 

^^ "■ 1^ f 

—m- i — *■ — \ ~~' ' ^ 

k+P P -^^ / 

(a) (b) (c) 



Figure 6: Two loop contribution of fermion pair interaction to wave function renormalization 
ofY, Y\ They contribute to I operators in the dilatation operator 

k + l 



va.^>v 




P k+l+p P 

(a) (b) (c) (d) 

Figure 7: Two loop contribution of vacuum polarization to wave function renormalization of 
7,7'*'. Both U(N) aiidU(N) gauge parts give additive contributions. 

diagrams involve parity-even vacuum polarization, as shown in Fig. U\ Their computations are 
summarized in Appendix B4. We also present the analysis of the one-loop vacuum polarizations 
in Appendix B5. 

Summing up all these wave function renormalization to 7, 7^, we find their contribution to 
the dilatation operator as 



Hz = l^ 



1 2 1\ /4 \ 3^ 2L 



12 + 3 + 3J + l3+^ 



Li 



= ^'E^I (4.15) 

In the first line, the first parenthesis is the contribution from gauge fields: diamagnetic inter- 
actions, paramagnetic interactions, and Chern-Simons interactions. The second parenthesis is 
the contribution from fermion fields. The last term is the contribution of vacuum polarization. 
Adding up all the contributions, 

^totai = HB+Hp + Hg^ug^+Hz (4.16) 

we get the result (|4.8I) . As claimed, this is precisely the alternating spin chain Hamiltonian 
we obtained from mixed set of Yang-Baxter equations. As such, we conclude that dilatation 
operator of iA^ = 6 superconformal Chern-Simons theory of ABJM is integrable at two loops. 



16 



We stress the importance of explicit and direct computation of the dilatation operator with- 
out a prior assumption relying on supersymmetry or integrability. It is satisfying that the result 
passes various compatibility tests. For instance, take chiral primary operators. These are subset 
of the single trace operators O where Y's and F^'s are totally symmetric and traceless under any 
contraction between 7' s and F^'s, and corresponds to massive Kaluza-Klein modes over CP^ in 
the Type IIA supergravity dual. Because of supersymmetry, their scaling dimension should be 
protected against radiative corrections. Indeed, acting on these operators, //total vanishes since 
contribution of terms involving K operator are null and contribution of P cancel against that of 
I. As a corollary, the fact that our result is consistent with expectation from supergravity dual 
implies that the dimension reduction method we adopted for computations are compatible not 
only with Slavnov-Taylor identities of the gauge symmetry but also with 9l =6 supersymmetry. 

5 The Shortest Chain and Wrapping Interactions 

In deriving the dilatation operator in the last section, we assumed that the gauge invariant oper- 
ator is infinitely long, L ^ oo. From planar diagrammatic s, we see easily that dilatation operator 
computed perturbatively up to the order 2i will give rise to a spin chain Hamiltonian whose 
range extends to (2£)-th order. Therefore, for operators of finite length, a new set of planar dia- 
grams which wraps around the operator will come in to contribute. These are so-called wrapping 
interactions, a feature discussed much in the context of integrability of four-dimensional 9{, =4 
super Yang-Mills theory [|l2l, [I44l-[|53]|. 

In 5V! = 6 superconformal Chern-Simons theory, the situation is more interesting. Since the 
dilatation operator at two loops ranges over three sites, spectrum of the shortest gauge invariant 
operator of length 2L = 2 will receive contributions from wrapping diagrams already at leading 
order! In this section, we like to identify these wrapping interactions for the shortest gauge 
invariant operators and discuss their implications. 

Let us denote basis of the shortest operators as 

1/1/2) = TvY^^yI = Oi,i, G 4®4 . (5.1) 

The 4 (g) 4 representation is decomposed irreducibly into the traceless part, 15, and the trace 
part, 1. The multiplet 15 is chiral primary operator, so their conformal dimension ought to be 
protected by supersymmetry. 

To check this, let us first identify the two-site dilatation operator that includes the wrapping 
interactions. At two-loop orders, the scalar sextet interaction does not contribute to length- 
2L — 2 operators since only four legs can be connected to the operators, leaving a tadpole that 
vanishes identically. Hence the dilatation operator consists of the two-site plus wave function 
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w 



Figure 8: Two loop wrapping interaction contribution to the shortest gauge invariant operators, 
(a) fermion field wrapping, (b) gauge field wrapping, (c) a new gauge triangle. 

renormalization parts plus wrapping contributions. 

From the computations of section 4, the original two-site contributions comprise of two-loop 
diagrams from gauge interactions and from Vp interactions. Their contributions are 



H2 



2 4 



x2 



K- 



l^w. 



(5.2) 



In the first line, the first term is the contribution of gauge interaction diagrams and the second 
term is the contribution of Vp interactions. We computed total energy, so multiplied the energy 
density by the spin chain volume 2L — 2. The one-site contribution arising from the wave 
function renormalization is 



Hi 



-l'^Ix2=-l'^I. 
4 2 



(5.3) 



Adding these two and acting on 15 in (15.11) . we see that anomalous dimension of the chiral 
primary operator is non- vanishing. If our regularization method of dimensional reduction plus 
infrared mass regularization were compatible with supersymmetry, there must be other contri- 
butions heretofore unaccounted that would cancel against the non-vanishing contribution (15.31) 
and protect the anomalous dimension of chiral primary operator from quantum corrections. 
These are precisely wrapping interactions. 

Indeed, for the shortest operators of L = 1 under consideration, there are three classes of 
nontrivial wrapping interactions. We now summarize their contribution and relegate details of 
Feynman diagram evaluation to Appendix C. 

There is the gauge field wrapping contribution with the diamagnetic interactions as in Fig. [8] 
(a). Its contribution is 



H. 



glw 



XH. 



There is also the fermion field wrapping contribution as in Fig.[8](b). Its contribution is 



H. 



2(K-I)X2 



(5.4) 



(5.5) 
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It is important to note that these two wrapping interactions utilizes simultaneously U(N) and 
U(N) interactions. Thus, this contribution arises not just by distinct topology of planar diagram 
but from very different interactions from the original, unwrapped two-site interactions. 

There is also a doubling-type wrapping contribution of using the same gauge group in- 
teractions. This happens only for the gauge interaction diagram contributing to K operator. 
Moreover, the contribution is doubled since there are two distinct ways of wrapping. This is 
best illustrated on a cylinder, from which we see that there are two different kinds of topology 
of wrapped Feynman diagrams. From appendix C, we identify this contribution as 

HgKw = -X^K. (5.6) 

Putting both the original and the wrapping diagram contributions together, the full Hamilto- 
nian of 2L = 2 operator is given by 

H2L=2 = 2l^K. (5.7) 

Notice that the part proportional to I operator is canceled between the original and the wrapping 
interaction contributions. One thus check that the chiral primary operators 15 indeed has a 
vanishing anomalous dimension since, by definition, it has no trace part and is annihilated by K 
operator. For the singlet 1, I*) = j\II), the anomalous dimension is 

H\s) = S'k}\s). (5.8) 

It is interesting to compare the above spectrum with spectrum of the naive Hamiltonian 
^naive, viz. the alternating spin chain Hamiltonian with periodic boundary condition and 2L = 2. 
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The latter is 

"naive ^ ^ 2^ 



2 r 1 1 



2X2k, (5.9) 



for 2L = 2. Acting on 15 and 1 states, we find that their anomalous dimension is and 4 ■ 2?i^, 
respectively. So far, we computed the spectrum of the shortest operators without a priori as- 
sumption of supersymmetry. As a consistency check, we now compare these spectra with their 
superpartners. Recall that length 2i operators with Dynkin labels {i — 2m, m + n,i — 2n) and 
length 2£ — 2 operators with Dynkin labels (£ — 2m, m -|- n — 2, £ — 2n) are superpartners each 
other. Here, we have the simplest situation: the L = 1 operator 1 of Dynkin labels (0, 0, 0) is the 
superpartner of L = 2 operator 20 of Dynkin labels (0, 2,0). Fortuitously, anomalous dimension 
of the latter was computed at two loops in [i33il to be 8?i^, and matches perfectly with our com- 
putation [^. Note that, at two loop order, the L = 2 operator 20 does not receive any wrapping 



'•^Here, we used the identities ( 13.21b . 

'^We thank Joe Minahan and Kostya Zarembo for useful correspondences on this issue. 
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interaction corrections. As such, we may consider agreement of the anomalous dimensions 
between the two superpartners as a nontrivial confirmation for the wrapping interactions we 
studied for the L = 1 operator 1. 

We should also note that the naive Hamiltonian is not the right dilatation operator for the 
shortest operators. Nevertheless, interestingly, the spectrum of naive Hamiltonian coincides 
with the spectrum extracted from the true two-site Hamiltonian. It would be very interesting to 
see whether this coincidence persists to higher orders in perturbation theory. 

6 Bethe Ansatz Diagonalization 

In section 3, we constructed transfer matrix. To obtain spectrum, we need to diagonalize the 
transfer matrices. Within algebraic Bethe ansatz, a fairly general result is known for a Lie 
(super)groups G [l30l[26l|29l. It suffices to adapt the results to the case that G =SU(4)H Dynkin 
diagram of SU(4), drawn horizontally, has three roots: left(l), middle(m), and right(r). The 
diagonalization is specified by the choice of Dynkin label (Ri.Rm.Rr) for the site representation 
R and total number of sites L^ that representation occupies. In the present case, we have placed 
4 and 4 representations at alternating lattices, so Ri = Rr = l,i?m = and L4 = L j = L. Each 
excitation is associated with three sets of Bethe ansatz rapidities {la.nih, rc)'s whose labels range 
over [IjA'z], [l,A'„j], [IjA'^r]^ respectively. It belongs to the SU(4) representation with the Dynkin 
labels {L-2Ni+Nm,Ni + Nr-2Nm,L-2Nr + Nm). Positivity of the Dynkin labels restricts 
range of the three Bethe ansatz rapidities accordingly. Then, choosing the highest- weight state: 

L _ 

|^+) = n®|l)2^-i|4)2^=|1414---) (6.1) 

as the ground-state, the eigenvalue of the transfer matrix Tq{—u) o is found to be 

A{u) = {u-l)\u-2)'^Yl ''l'^''^\ +u\u-l)'^Yl ''~J'''~_l (6.2) 

a=l U — lla— 2 c=l ^ ~ ^^c ~ 2 

+ ..^(„-2)^[n"-''°-fn"-''""--°+n"-''"^-^n"-''--ii. 

\-^J[ u - ila - J hJ[ u-inib-l ^J[ u- lint- I ^Ji u-irc-^^ 

We have chosen the Bethe rapidities symmetric between the three roots. Keeping the highest 
weight state the same 1 1414 ■ ■ •), we also find that diagonalization of the second transfer matrix 
Tq{—v) proceeds much the same way as that of 7b (—v) except that we interchange role of the 



'^For SU(3) alternating spin chain, this was done explicitly in |28|. 

'^For later convenience, we choose to diagonalize Tq for opposite sign of the spectral parameter u. 
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left and the right SU(4) roots: 

_ Ni 

A(v) = v\v-\fY[- 



iln 



Nr 



Ni 



il -^ 



+(v-i)v-2)^n 



■ire +2 



+ v\v-i)^ n 



V - ila 



T 

2 b 



c=\ 

N„, 



n 



i V - //a - f ^Ji V - imt - 



T\n 



V — zr.. 



V — inii, 



(6.3) 



iVr 



^ V — inib 



-n ^ 

1 -•■-•■ V — ;V — 1 



3, 

2 



Mutually commuting conserved charges are then constructed by expanding these eigenvalues 
around m, v = 0. The first two charges are the total momentum and the total energy: 

1 



^total 



logA(M)+logA(M) 



M=0 






rb + '/2 
n - '/2 



(6.4) 



£^total 



X^ 



X^ 



— (logA(M)+logA(M)) 
dw 



u=Q 



1 



A'/ ;^ A'r 



(6.5) 



Here, we chose fundamental domain of the momentum to [0, 2%) and scaled the total energy by 
y? in accordance to the relation we fixed between Hamiltonian derived from Yang-Baxter equa- 
tion and from superconformal Chern-Simons theory. Likewise, we can deduce higher conserved 
charges from higher moments of the transfer matrices. 

The Bethe equations that results from the above transfer matrix eigenvalues A, A are 

L 



/ -- 

fl 2 

'a ~r 9 



T-x 'a l^b ' "TT" '•a '"c T^ 2 



n 

b=l{b^a 



nin 



■nib + i 



nia-nib 



ra + i 



n 

b=l{by^a) 






N,n 

n 



m. 



f —i Ny 

/"a - mc + 2 



m« 



■^rf 



nin 



rd + h 



-A 'a 



-nir — 



(6.6) 



It is straightforward to check that these same set of Bethe ansatz equations remove potential 
simple pole terms for both A and A simultaneously. 

From the integrability perspectives, (2 -|- 1) -dimensional superconformal Chem-Simons the- 
ory is quite different from (3 -|- 1) -dimensional super Yang-Mills theory. The most distinct fea- 
ture is that the spin chain associated with dilatation operator is not homogeneous but alternating. 
It calls for better understanding to questions that arise in comparison with fA^ = 4 super Yang- 
Mills counterpart. We shall now study spectrum of the Bethe ansatz equations for a few simpler 
situations and gather features concerning excitations of the alternating spin chain system. 
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First, consider the special class of A^„, = for arbitrary L>2. The first and the third Bethe 
ansatz equations decouple, and each equation becomes the same as the Bethe ansatz equation 
of the well-known SU(2) XXX i spin chain. Thus, if one can identify the first set with SU(2) of 
4 side, then the third equation corresponds to SU(2) of 4. We then have two decoupled sets of 
the solution including towers of bound states, and they are exactly the same as the XXX i spin 
chain. 

Now let us consider the case Ni = Nm =Nr=\ case for a general L>2. The Bethe ansatz 
equations are reduced to 

L 

/ — 4 \ l — m + i 



1+^ I l — m — ^ 

m = -{l + r) 

r—i\ r — m+ 1 

— T = T- (6-7) 

r+jj r-m-\ 

In terms of the individual momentum variables, after using the second equation, the combina- 
tion of the first and the third equations becomes 

This is solved by 

27r?2 
P = Pi+Pr = —. (n = 0,l,2,---,L-l). (6.9) 

First, consider m = case. In this case, total momentum P = 0. For the relative momentum 
q = {pi — pr), we also have 



This is solved by 



ML+i)/2^^_ (6.10) 



^ = ^^. (6.11) 

2 L+1 ^ ^ 

For this case, the energy (|6.5I) is given by 

E = Sl?sm^-. (6.12) 

This is the simplest example of two-particle excitations where 4 and 4 excitations are correlated. 
The total momentum is zero, while total energy depends on relative momentum. 
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Consider next general m. From the ratio between the first and the third equations, we obtain 

;9r l — r — i .^ ,^. 

l — r-^i 

where we used the second Bethe ansatz equation for simplification. Total momentum P is 
nonzero. Furthermore, expressing this equation in terms of P and q, we find the relations: 



P sinf(L + 2) cosf(L + 2) 

— =^ or 

2 siniL cosIL 



cos- = ^^- — or ^^- — . (6.14) 

^ ""l^^ cos^/ 



If q is real, viz. two real Bethe roots, the relation shows that relative momentum q is correlated 
with total momentum P. That is, even though there are two excitations associated with 4 and 4 
chains, their motion exhibits mutual correlation. If q were imaginary, viz. a Bethe string, the 
relation shows that total momentum ought to be purely imaginary. This show that there cannot 
arise any bound-state between 4 and 4 spins. 

We can also comment on thermodynamic limit in which densities of the Bethe roots are kept 
finite. By taking L-^ oo limit of the Bethe ansatz equations and taking the so-called "no hole" 
excitation condition, we obtain relations among the three Bethe root densities p/ (x) , Pm(-f ) , Pr{x) . 
From the first and the third Bethe ansatz equations, after Fourier transform, we find 

pi{k)=pr{k). (6.15) 

This has a simple interpretation: because the alternating spin chain is manifestly charge-conjugation 
invariant, excitations ought to be so as well. Moreover, from the second Bethe ansatz equation, 
we obtain 



P™(fc).-W/2 = l[p,(fc)+p,(fc) 

It immediately follows from these two equations that the mean value of root densities 



(6.16) 



We conclude that all three Bethe root densities are equal, and hence 4's and 4's are equally 
populated and balanced each other for the minimum energy configuration. 

Furthermore, analysis of the shortest operator suggests that excitation in superconformal 
Chem-Simons spin chain is different from excitation in iV; = 4 super Yang-Mills spin chain. 
In the latter, the vacuum is ferromagnetic and excitations break 80/^(6) to [SU(2)]^. The latter 
is the symmetry group of dilute, finite-energy excitations. In the present case, analysis of the 
previous section seems to indicate that excitation is organized by the full SU(4), not by any 
subgroup of it. This is because the finite energy excitation is a singlet of SU(4), not of any 
subgroup of it. Lastly, in this system, excitations with Nm = comprises of two decoupled 
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XXX 1 spin chains with its own ferromagnetic vacuum, respectively. Though this is certainly a 

1 

closed subsector, general excitations in the full system looks quite different, as is seen above in 
the simple situation of Nm = 1 . 

Following the general prescription ||29l and paving the parallels to what was done in the con- 
text of fA(; = 4 super Yang-Mills theory jSl, extending the SU(4) spin chain to the OSp(6|2, 2; M) 
superspin chain and writing down Bethe ansatz equations are immediate and straightforward. 
This was done already in [|33l . More recently, spectrum in the Penrose limit [[54ll . various 
SU(2|2) closed subsectors 11311 . all loop Bethe ansatz equations ll55l . and finite-size effects [f56l 
were studied. With these developments, it would be interesting to explore precision tests for the 
new correspondence proposed by ABJM. 
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A Notation, Convention and Feynman Rules 
A.l Notation and Convention 



^^'^ metric: 



g,„„ = diag(-,-F,-F) with m,n = 0,1,2. 

e012 „ ,1 

e = — £012 = +1 



(A.l) 



L^'^ Majorana spinor and Dirac matrices: 

\|; = two-component Majorana spinor 
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x^« = £«P^p, ^„ = £^p^P where £«P = -£„p = /o^ 

7SP = (/o2,o3,oi), (y")„p = (-l,oi,-o3) Obeying ^f = r"-e""%-(A.2) 

A.2 ABJM Theory 

• Gauge and global symmetries: 



gauge symmetry : U(N) (g) U(N) 
global symmetry : SU(4) 



(A.3) 



We denote trace over U(N) and U(N) as Tr and Tr, respectively. 

• On-shell fields are gauge fields, complexified Hermitian scalars and Majorana spinors 
(7=1,2,3,4): 



A«: Adj (U(N)); A„ : Adj U(N) 

Y^ = {X^ + iX^X^ + iX^.X^ - iX\x'^ - iX^) : 
y/ = {X^-iX^,X^-iX^,X^ + iX\X^ + iX^) : 
^i = (¥^ + a^ -¥^ - a\ ¥4 - ^4, -V3 + as) : 



(N,N;4) 
(N,N;4) 
(N,N;4) 
(N,N;4) 



(A.4) 



action: 
I = 



1.2 



k ( 2i 

— £ Ir I AffidfiAp -\- —AffiAfiA 



k 

4%' 



£ Ir I Afjid„Ap + —AffiAnAp 1 



-Vf-Vb 
Here, covariant derivatives are defined as 

Lfml ^ 'Jml "T 'Amy ^y ^m 7 L/ynli ^ '^mli "T zA^yy zy^ A^ 

and similarly for fermions ^/, ^^^. Potential terms are 



(A.5) 



(A.6) 



27rz— 
Vf = ^Tr 



y^ty/vj/tyvj/^ _ 2y/y^vj/t/ij/^ + £^-^^^y/^yy^^i] 



2%i 



Tj.[y/y^tvj/^vj,t/ _ 2y^y/^/^t/ ^ £^^^^y^vi/t7y^vi/tL 



(A.7) 
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and 



- - -m- 



2 

y^ty/yty/f y^y/ ^ y^ty/yty/y^y/T 



^4y^ty7y^y/yty/^ - 6Yj'Y'y}Y^YIY-' 



(A.8) 



At quantum level, since the Chern-Simons term shifts by integer multiple of Stt^, not only A'^ 
but also k should be integrally quantized. To suppress the cluttering 2n factors, we also use the 
notation K = J^. At large A^, we expand the theory and physical observables in double series of 

1 A^ A^ 

by treating them as continuous perturbation parameters. 

A.3 Feynman Rules 

• We adopt Lorentzian Feynman rules and manipulate all Dirac matrices and Emnp tensor ex- 
pressions to scalar integrals. For actual evaluation of these integrals, we shall go the Euclidean 
space integral by the Wick rotation, which corresponds to x^ -^ —ix. In the momentum space, 
this means we change the contour of po to the imaginary axis following the standard Wick ro- 
tation. Then in terms of integration measure, we simply replace d^^k — > id^'^kE together with 
p^ -^ +p\- The procedure is known to obey Slavnov-Taylor identity, at least to two loop order. 

• We choose covariant gauge fixing condition for both gauge groups: 

d'"Am = and a'"A„ = (A. 10) 

and work in Feynman gauge by setting the gauge parameter ^ to unity. Accordingly , we 
introduce a pair of Faddeev-Popov ghosts c,c and their conjugates, and add to / the ghosts 
action: 



^shost 

^ '^2,1 



Tra'"c*D„c + Tra'"c*D„c 
Here, D^c = dmC + i[Am, c] and D^c = dmC + i[Am,c]. 



(A.ll) 



Propagators in U(N)xU(N) matrix notation: 



r 



gauge propagator : A„„ [p) = —I — — - 

k p^ — i£. 



scalar propagator : D/{p) = 5; 



— z 



p^ — ie 



fermion propagator : S' j{p) = ^j —5 — 



p^ — ie 

ghost propagator : K(p) = -^ (A. 12) 

p^ — z£ 
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• Interaction vertices are obtained by multiplying / = ^/^ to nonlinear terms of the Lagrangian 
density. Note that the paramagnetic coupling of gauge fields to scalar fields has the invariance 
property under simultaneous exchange between Am,!'^ and A^^Yj . 

B Two-Loop Computations 
B.l Two-loop integrals 

We first tabulate various Feynman integrals that appear recurrently among two-loop diagrams. 
They are all evaluated straightforwardly by Feynman parametrization 



A«5* r{a)r{b)J J ' \Ax + ByY+''' 

We use dimensional regularization by shifting the spacetime dimension to d = 2(D = 3 — e. The 
ultraviolet divergence shows up as a simple pole l/e. It is related to the momentum space cutoff 
A as 

-:=21ogA. (B.2) 

In the following, we collect factors arising from propagators in parenthesis and those from 
vertices in square bracket. We have the following integrals: 

^2co^ ^2co£ J J ^ 



h = 



(27l)2«(27l)2»(;t + £)2/t2£2 

1 . r (f^i 1 r (f'^k 1 

ax 



J (27r)2to£2y (27r)2co[yt2 + 2;c/t-£+;c£2]2 

1 /•! 1 r (f^^i 1 

ax- 



1 1 1 

The integral that appears in fermion and gauge boson exchange diagrams is: 

^ f d'"fe d2">£ 1 2{k + i).i 
^ J (27r)2»(27r)2«(/t2)2(yt + £)2£2- '^ • ^ 

We perform the £ integral first after using the Feynman reparametrization: 

J _ [\ [ ^^"^k 1 /■ d2"£ 2£-(£ + /t) 



J {2n)^^{k^)^J {lnf^[fi + 2xk-l+xk^Y 

1 /•! v^ r d2"/t 1 

--— ax ' 



SnJo v^r^y (27r)2»fc3 

1 1 1 
8 47r2i' 
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(B.5) 



where for the second equality, we have used the integral, 



(f'^i 






1 



x^k^knr{l/2) gmn r(-l/2) 



+ 



(27r)2»[£2 + 2;c£-fc + fc2]2 (471)3/2 L [4 l-;c)/t2] 1/2 2 [x{l - x)k^]-^/^ \ 

(f^i tyn _ 1 xk^T{\IT) 

(271)20) [£2 + 2;c£-fc + fc2]2 - -(4^)3/2 [_^(1 -;c);t2]l/2 " 



(B.6) 



If one exchanges the order of integrations, there may appear an infrared singularity. How- 
ever by introducing infrared regulator mass m, one may get the same result in the limit (0^3/2 
and m ^ 0. 

In the gauge boson exchange diagram, the following integral appears: 

d2"/t /■ d2"£ 



h 



1 



1 1 



(27r)2« y (27r)2» (/t + £)2 {k^Y (£2)2 
We evaluated them as follows: 



{{k-tY-kH'-\=h^-hs 



Ha ^ 



d^^k k^kn f\^0) f d^'^^ il-x)emL 



(27r)2» (/t2)2 70 r{2)J {2nY^[i^ + 2xi-k + k^Y 

d^^k k k r^ 1 

dx(l — X 



(27r)2» (/t2)- JO 

1 1 1 



(471)3/2 L 



■^t^K ^P/^' 



+ 



r(i/2) 



[4l-;c)yt2]3/2 2 [x(l-;c)/t2]i/2j 



-^3,^ 



16 47^2 £■ 

d2«) 1 .1 



(27r)2»/t2yo 
r(l/2) rl 



da 



dx- 



j2co/ 



1 



(27r)2»[£2 + 2;c£-/t + ;c/t2]2 
1 r d^'^k 1 



Hence, 



(47r)3/2r(2) Jo X{1-X)J (271)2" (/t2)3/2 
1 1 1 

847r2£ 



_ 1 1\ 1 1 _ 111 
^^~ 'T6~8v'47r2£~~T6 47r2£' 



(B.7) 



(B.8) 



B.2 Contribution from Sextet Scalar Potential 

The Lagrangian contains sextet scalar interaction — Vscaiar- Three of the scalar fields couple to 



and the rest three to 0^. With U(N) and U(N) index loops, combinatorial factors are given 
by 



3-A^^ 



IT®' 



2r - 4Pi3 ® I2 - Ki2 ® I3 - 1 ® K23 ® I2 + 2Ki3 ® K12 + 2Ki2 ® K13 



(B.9) 
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There are three scalar propagators and one interaction vertex, contributing factors 

1 Q O 471^ T 

,(/)3[/].A^2_ ^2 (B.IO) 



3k2^"' '-' " 3 

The remaining 2-loop integral is given by /i. Summing over all contributions, the scalar sextet 
interaction gives rise to 2-loop dilatation operator 



^ - ' '- - - " - ~ (B.ll) 



^ £=1 



B.3 Contribution from two-site Interactions 

In this appendix we shall present the full detailed computation of the two site interactions. First 
let us compute the Yukawa two-site interactions. The nonvanishing Yukawa interaction leads to 
only a K-type interaction. The relevant Feynman diagram is depicted in Fig.|2]D. 

With two Yukawa interaction components and one U(N) and one U(N) color traces, combi- 
natorial factors are gathered as 

— ■2-N^=N^. (B.12) 

2! 

There are four propagators and two vertices. This yields numerical factors 

i-i)Hin?{±^) •(-)fd (B.13) 

where the subscript ( )fd signifies the Fermi-Dirac statistics minus sign. The loop integral is 
given by 

.2 f (f'^k d2"£ 1 // Jt+jt 



where the f factor comes from the analytic continuation of the integration measure. 

After taking the gamma matrix trace try^y = 2g""\ this integral equals to —I2 in (IB. 51) . 
Hence putting everything together, one has 

^2^ ^K (B.15) 

2£ 

for the Yukawa two-site interactions. The contribution to the operator renormalization is nega- 
tive of this: Therefore, the Yukawa contribution is 

2L 

//f = ?i'£K^,£+i (B.16) 
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We now evaluate the gauge two-site interactions. 

The gauge boson interactions contribute both IK and I type diagrams to the dilatation op- 
erator. Let us begin with IK type contribution. The relevant diagram is in Fig. [2t. It has 
combinatorial factors 

— ■2-N^=N^. (B.17) 

There are three boson propagators, two gauge propagators and one seagull interaction vertex. 
So, numerical factors are given by 



where the last factor accounts for the (±) relative sign of U(N) and U(N) Chern-Simons term. 
It is important to note that the gauge field propagator in momentum space has no / = -\/— T. The 
loop integral reads 

where again the i^ factor comes from the Euclidean rotation. Using the identity g"^Zmnp^qrs = 
—(gnrgps—gnsgpr), wc find that the integral is the same as 4/3. 
Hence putting everything together, one has 

^ '-K (B.20) 

2 2£ ^ ^ 

for the gauge two site IK contributions and, for the operator renormalization, 

T, 2 1 

K. (B.21) 

2 2£ ^ ^ 

There are also contributions to I from ^-channel exchange of diamagnetic gauge boson inter- 
action. The corresponding Feynman diagram is depicted in Fig [2^. There are two scalar prop- 
agators, two gauge boson propagators and two diamagnetic vertices. Note again, for Chern- 
Simons theory, gauge boson propagator has no i in momentum space. So, the combinatorial 
factor is 

1^2-i-i)^-[i]^-N^-(^y = {ATi^)X\ (B.22) 

The loop integral reads 



•i 



(27r)2« (27r)2to (/t2)2 (k + i)^ e 
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(B.23) 



Using the identity £™"'^£„„j* = —Ig"^, we find that this integral is the same as —h- There 
are identical contributions from each letter (with alternating U(N) and U(N) gauge boson ex- 
changes), we find the contribution as 

4<^.. (B.24) 

The corresponding operator renormalozation contribution is 

X^ 1 

I. (B.25) 

Alt ^ ' 

Hence there are two gauge two-site contributions. Using l/(2£) = InA, the gauge two-site 
contributions to the anomalous dimension are summarized as 



"gauge — 2^ 



2L 

l^. (B.26) 



il-lK,,^+i 



B.4 Contributions of Wave Function Renormalization 

The first one involves diamagnetic gauge interactions. The relevant Feynman diagrams are in 
Fig. [3l As scalar fields are bi-fundamentals, there are processes involving U(N) gauge boson 
pair, U(N) gauge boson pair, and one U(N) gauge boson and one U{N) gauge boson pair, which 
are respectively corresponding to Fig. [3^, Fig. [3h and Fig. [Sj:. Taking account of opposite rel- 
ative sign between gauge boson propagators for U(N) and U(N) and of different combinatorial 
weight of diamagnetic coupling terms, the numerical factor reads 



-2 • i-i) . [i]^. [(-)2 . (+)2 + (+)2 . (_)2 + (+)(_) . (_2)2]A^2 (^_) = -2i(4n')l^ (B.27) 
Denote momentum of the external scalar field as p™. Then, loop integral reads 



(f^k d^'^i 1 e"""A:a£,„/ 



' J (27r)2«(27r)2»(A: + £ + p)2 k^ f ' ^ ^ ^ 



(27r)2« (271)2" (A: + £ + p)2 F 
For the evaluation of this integral, let us introduce 

^20)^ ^2C0£ J 2/t 



hip) 



i2n)^^{2n)^^{k + i + p)^k^fi 



d2"£ 1 /■! r d^'^k 

dx 



{2n)^^i^Jo J (27r)2»[(/t + x(£ + p))2+4l-;c)(£ + ;?)2]2 
^ f'dx r^ f d'"^li£M (B29) 
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The x-integral is finite and equals to n/2. The remaining ^-integral can be performed by 
applying Feynman's parametrization. In dimensional regularization, we have 

ir(3/2) /-i 1 /■ d2«^ l-p 



:3/2) /•! J_ r _6^ 
[1/2) Jo \/yJ (2n 



8 r( 1/2) 70 \/yJ (27r)2»[£2 + 2x£-p+x;?2]3/2 



1 /-i dy 



r(£) yp 



2 



16 70 v^L {ATirr{3/2) {y{\-y)p^Y 
Taking e = 3/2 — co ^ 0, this integral equals to 



(B.30) 



Putting together, we thus find that these diagrams contribute to the wave function renormaliza- 
tion as 

-Ix2l(/p2) (B.32) 

Consider next two diagrams involving four paramagnetic couplings. Planar diagrams in- 
volve two vertices from U(N) and two from U(N), as shown in Fig. HI Taking care of opposite 
relative sign of gauge boson propagators between U(N) and U{N) and that there are three inter- 
nal scalar propagators, we have combinatorial factors 

^22 . (+)(-) . (-0^[/]'(^) V(2) = -2i{4n')l\ (B.33) 

where we put an additional factor two because there are two such diagrams. With external 
momentum /)"\ the loop integral read 

^2co^ d2"£ £'""^{i + 2p)^{p + 2k + 2i)niqe'''''{2i + k + 2p)a{k + 2p)bkc 



i- 



(B.34) 



(271)2'^ (27r)2» (^k + e + py{i + p)^{k + p)^kH^ 

This integral can be integrated without further assumption but we note that the numerator of the 
integrand is already quadratic in pm- Using the isotropy of the system, we replace 

p2 

PaPb -^ -:rgah (B.35) 

and then set p to zero in the remaining integral. One may show that the results from the both 
methods agree precisely with each other. 
Thus the integral becomes 



Putting all the factor together, one has 



-^?^'^0'/) (B.37) 
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There are two diagrams involving Chern-Simons cubic coupling. The contributions of U(N) 
and U(N) are added up with an equal weight. The Feynman diagrams are in Fig. [51 
The relevant combinatorics is 

3!-3, ,,2r..l4^.2/(±)K^/(±l)^^o^ o.v.^2^.2 



3! 



i-i)^[i]^N'^^) (m) (2) = -2i{4n')l' (B.38) 



where the last factor two takes care of the U(N) contribution. The loop integral becomes 



•2 



d^^k d2"£ e"''l%nar{2p + krk'-£nbsi2p-i + k)''i-k-i)%cti2p 



\bnt 



(27r)2» (271)2^0 (p _ £)2(^ + £)2(^ + pY^2(2 

Using the rule of 



(B.39) 



the integral becomes 



Using /3, one finds 



^ab 2 
PaPb -^ -^P 



.28^2 f (f^k (f'^i -k^f+{k-if 

[2%)^^ [2%)^^ {k+iY{k'^Y{eY 



p2 1 1 



(B.40) 



6 47r2 £ ■ 
Therefore, the whole contribution combining the combinatorics becomes 



(B.41) 



-ix2l/p2 (B42) 

There are also diagrams involving paramagnetic and diamagnetic couplings. Their net com- 
binatorial factor is nonzero, but the loop integral vanishes identically. 

Let us now turn to the Yukawa contributions. First consider the Feynman diagrams in Fig. [6^ 
and Fig. [6J5. Within the planar diagram, both fermion can be joined either U(N) side (Fig. [6^) 
and U(N) side (Fig. [6b). The joining using the first two terms has a factor 4 from the SU(4) 
index contraction. Then the cross terms between the first two and the second two terms in total 
have a factor —4. Hence one can check that this cross contributions cancel precisely the those 
from the first two. 



By combining the second two of the Yukawa potential, for the U(N) and U(N) side, we have 
combinational factors 

i2-(-/-/2)[/]2.(^|jV(-)FDx8 = -32/(47r2)X2. (B.43) 

where the extra factor eight comes from one contraction of SU(4) index and the doubling by 

U(N) and U(N). 
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Then the remaining integral has the expression, 



which is same as Iq- Therefore the whole contribution is 

-^X2l(/p2). (B.45) 

For the wave function renormalization, the third two of Yukawa potential also contribute. 
The diagram is in Fig. (6]:. It has combinatoric factors, 

(2!)2 . (-/) . (/)2[/]2 . (1) (^) A^'(-)fd X (-6) = -24i{4ii')X' , (B.46) 

where (2!)^ is the usual symmetry factor of the Feynman diagram. The last factor (—6) comes 
from the following SU{4) index contraction 

eiABce-'^''^ = -6 8i (B.47) 

where / is for the incoming and the J for the outgoing scalar SU{4) indices. 
Then the remaining integral takes precisely the same from: 

2 f d^'^k d^'^l I trjtk 

i^ / .._^...w._^...w . , .^.,.^; (B.48) 



(27r)2" (27r)2» (p + k-i^k^i^ 
which is again the same as Ig- Therefore the whole contribution is 

-l^-iip^). (B.49) 

Finally, there are the vacuum polarization contributions of the gauge loop. The relevant 
diagrans are depicted in Fig. U\ 

As we shall explain in the following appendix, the self energy correction for both A and A 
gauge fields is given by 



iUatik) = 8/ 



'^a'^b §ab'^ 



(B.50) 



I6k 

where the factor eight comes from the four complex scalars and fermions with an equal weight. 
For the relevant diagram of Fig. U\ the combinatorics factor reads 

|-(-0['f • Q Vx (2) =2/(47r2)X2, (B.51) 

where the last factor two comes from the doubling by replacing A gauge by the A gauge field. 
The remaining Feynman integrals takes the from, 

d^^k e^'^"{2p + k)^{-k)„eb'J{k + 2p),kjiUab{k) 



d^'^k k^p^-{k-pf 



= 2/^ 

J (271 



(271)2'^ {k + pfk^ 

34 



(B.52) 



where we have a single i produced by the Euclidean rotation. 
By the dimensional regularization, this leads to 

1 1 



""We- <«■"' 



Hence the total contribution reads 



All the remaining diagrams, one may prove that their contribution is identically zero after 
the dimensional regularization. 

Finally we add up all the above contributions to the wave function renormalization and find 
that 

-l^'l{-){-ip^). (B.55) 

Since the counter term is a negative of this, the two-loop scalar wave function renormaliza- 
tion becomes 

Z,^l--l^- = l--l\2\nA). (B.56) 

I 
In order to get the operator renormalization factor, one has to take Z/ out for each site, which 

corresponds to adding — ^ of (|B.56I) to the interaction part of renormalization. The final contri- 
bution to the anomalous dimension is 



2L 



12 33/V3 / 3 



2L 



3, 



»^ = ^T (l^ + ^^)+U+0-'h = ^\E?- <B.57) 



For the gauge two-loop contributions 1/12, 2/3, 1/3 including the gauge self-energy cor- 
rection contribution 8/3, the two-loop Feynman diagram computation is carried out in Ref. [l57l 
for the t/(l) case. One can check the precise agreement after taking care of the planarity fac- 
tor and the number of matter degrees. Furthermore, Ref. [|58l deals with the two-loop Yukawa 
contribution to the scalar wave function renormalization for again t/(l). This result is again 
matching with ours if one takes care of the planarity and the number of fermions. 

B.5 One loop self energy correction to the gauge field 



The self-energy correction enters in the same form for the U(N) and the U(N) gauge fields. 
Therefore we focus on the correction to A gauge field only. At the one-loop level, the boson, 
the fermion, the gauge and the ghost loops may in general contribute to the gauge self-energy 
correction. In this appendix, we identify these self-energy contributions. 
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We begin with the scalar loop contribution. It is the sub-diagram of Fig.|7^. The momentum 
k plays the role of the external momentum. The self energy contribution reads 



iK,{k) = {iy[in4)i 



d^"£ {2i + k)ai2i + k)h 



(In 



|2<B 



(yt + £)2£2 



(B.58) 



where the extra factor 4 comes from the fact that 4 complex scalars are coupled to the gauge 
field. Using the dimensional regularization, one obtains 



\6k 



iKbik) = (4)1 
Similarly, for the fermion loop, the self -energy contribution becomes 

m/m-m2[/i2mr i if ^'"^ TrY,(/+^)Y,/ 



(B.59) 



(B.60) 



where again the extra factor four comes from the fact that there are 4 complex fundamental 
fermions. Using the y matrix identity and the dimensional regularization, the contribution be- 
comes 



r/ 



i^ibik) = (4)/ 



\6k 



(B.61) 



Hence, each complex matter contributes by the same weight and sign. 

One can continue the dimensions 2(0 to four and obtain the vacuum polarization in four- 
dimensional Yang-Mills theories. The integration leads to the logarithmic divergence in this 
case contributing positively to the |3-function of the Yang-Mills coupling. Again, boson and 
fermion contributions add up. 

For the gluon self-energy contribution, we have 



i^ibik) = {3)-{3)[f 



IK 



1 

KJ 



T T /" d^' 



d2"£ e'"^"eJ"%^,enjr{i + k)'ii 



^2co 



(yt + £)2£2 






d2»£ (£ + k)aib+{i + k)bia 



^2co 



{k + £)^fi 



(B.62) 



It becomes 



ini,{k) = -i 



kgkb+ga bk 
32k 



2-, 



(B.63) 



which alone does not respect the gauge invariance. However, there exists also the ghost loop 
contribution. 



<ik)={in?i-)il^ 



d2«£ (i + k)aib+ii + k)bia 



(27r)2» 
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{k + * 



(B.64) 



where we put the extra (— ) sign due to the ghost statistics. Therefore, the ghost contribution 
cancels out precisely the gauge loop contribution, reproducing the well-established result [431 ■ 
Again, analytically continuing to four dimensions, the integral expression for the gauge part 
changes while the ghost integral remains intact. With Yang-Mills couplings, both contributions 
no longer cancel each other but contribute negatively to the [3-function. 

C Wrapping Interactions for the Two-Sites 

As in Fig. [HI there are occuring three kinds of wrapping interactions. First is the gauge interac- 
tions of two diamagnetic couplings in Fig. [8l). It is an I type interaction and happens, not for 
each site, but just once. 

The combinatorial factor is 

^2. {-if. [2i]'. {+){-). N^.{^y = -{4)4ii\\ (C.l) 

where {+){-) accounts for the the relative U{N) and U{N) Chem-Simons term and the factor 
two in the vertices takes care of the diamagnetic interaction. 
The loop integral 



(C.3) 



(27r)2« (27r)2» (A:2)2 (k + i)^ £2 
is the same as (IB. 231) . So, the loop integral is evaluated as 

1 1 1 

847r2£' 
Putting things together, we find the whole contribution as 

,2(-l)^ 



X^^ — ^I. (C.4) 

2e 

The corresponding operator renormalization contribution is 

X^^I. (C.5) 

The second is for the K type gauge wrapping, whose Feynman diagram is depicted in Fig. [8]:. 
It is doubling of the K-type interaction discussed for the general two-site gauge interactions. 
This doubling occurs due to the fact that, on the cylinder, one may have two different topology 
of the diagrams. Namely the diamagnetic interaction of the same gauge group may happen 
either one side or the other side, which is not possible for the infinite chains. From the previous 
result, the corresponding extra operator renormalization contribution is 

-^2^K, (C.6) 
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where we take into account of the fact that the doubling occurs both for the U{N) and U(N). 

There is an additional wrapping interaction coming from the third two terms in the Yukawa 
potential. The Feynman diagram is in Fig. [H In order to have proper contractions, one has to 
join operator site one (7^1) to the site two (Y'^^) whereas the operator site two Yj to YJ . The 
corresponding £ tensors in the Yukawa interaction produce 

ew2fie'2«^i^ = 2(I-K). (C.7) 

The combinatorial factors are gathered as 

^' ■ 2 . 2(-/)2(/)2[,f . N^ . (L) fzl] . (_)pj3 = 4(4;r2) ■ }} . (C.8) 



2! ^ ' ^ ' ^ J y^j y ^ 

The loop integral is given by 



(CIO) 



which is the same as (IB.14|) . By the loop integration, one gets 

1 1 1 

847r2i- 

Hence putting everything together, one has 

n}^ — ^(K-i) (c.ii) 

for the Yukawa wrapping interactions. Therefore, the Yukawa contribution to the operator renor- 
malization is 

2^2— (K-I). (C.12) 

Adding up the gauge and Yukawa contributions, the wrapping interaction contribution to the 
two-site Hamiltonian is 

//^,3p = I-K + 2(K-I) = -I + K. (C.13) 
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